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The runoff characteristics and harmonic analysis of the soil moisture 
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Abstract: Robinia pseudoacacia stands act as a typical ecological protection forest in hilly semi-arid area of China. Two fields 
of surface runoff were separately set up in R. psoudoacacia stand and its clearcut area in the western Liaoning Province 
(18°50’-122°25’ E, 40°24’-42°34’ N) for measuring the characteristics of runoff and sediment as well as soil moisture dynamics. 
Contractive analysis of the two land types showed that there existed a significant difference in volumes of runoff and sediment 
between the sites of R. pseudoacacia stand and its clearcut area. The runoff volume and sediment volume in clearcut area were 
much bigger than those in R. pseudoacacia stand, with an increase amount of 40%-177% for runoff and 180%-400% for sedi¬ 
ment. Hydrograph of surface runoff of typical rainfall showed that the peak value of runoff in R. pseudoacacia stand was de¬ 
creased by 1.0-2.5x10' 3 m 3 -s'' compared with that in its clearcut area, and the occurring time of peak value of runoff in R. 
pseudoacacia stand was 10-20 min later than that in its clearcut area. Harmonic analysis of soil moisture dynamics indicated that 
the soil moisture in R. pseudoacacia stand was 2.3 % higher than that in clearcut area, and the soil moisture both in R. pseu¬ 
doacacia stand and its clearcut area could be divided into dry season and humid season and varied periodically with annual 
rainfall precipitation. It was concluded that R. pseudoacacia stand plays a very important role in storing water, increasing soil 
moisture, and reducing surface runoff and soil erosion. 
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Introduction 

Robinia pseudoacacia is an important tree species to 
improving ecological environment, controlling soil erosion, 
adjusting river water table, and providing the fuel in Chi¬ 
nese “three-North” area (Wang 2002). Many researches 
have been carried out on hydrographic and ecological 
benefit of R. pseudoacacia stand in Loess Plateau (Wang 
1994; Zheng 1995; Wang 1991), and most of the re¬ 
searches focused on improving soil physical properties, 
increasing soil infiltration, decreasing surface runoff and 
soil loss, etc.. However, few researches were carried out on 
the water storage, water and soil conservation, and soil 
moisture dynamics in the artificial protection forest of R. 
pseudoacacia in hilly semi-arid area in China. 

The hilly semi-arid area is one of the important land types 
of the arid area in northern China, with an area about 
1.1x10 s km 2 . It is not only one of the key production bases 
of foodstuff, fruit and stock farming, but also the northern 
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doorway and ecological wall of North China and its south 
region (Liu 2000). In this region, there are many artificial 
protection forests of R. pseudoacacia, with the total area of 
1.4x10 4 km 2 . Therefore it is necessary to study quantita¬ 
tively the characteristics of runoff and sediment, soil mois¬ 
ture dynamics in R. pseudoacacia stand, and its clearcut 
area through a fixed observation. The aims of this study is 
to provide scientific data for regeneration, pattern ar¬ 
rangement, and comprehensive benefits assessment of 
artificial protection forest of R. pseudoacacia, and enrich 
and improve quantitative assessment theory about hydro- 
graphic and ecological benefit of protection forest in 
“three-North” region (Wang 1998). 

Study area 

The experiment was conducted at a typical hilly semi-arid 
area in the western Liaoning Province (18°50’-122°25' E, 
40°24'-42°34’ N). The experimental area has semi-arid 
seasonal climate, with annual mean precipitation of 
450-580 mm, annual mean evaporation of 2 088.5 mm, and 
annual mean temperature of about 8.5°C. The frost-free 
period is 142 days (Yan 1988). Vegetation of this area be¬ 
longs to north China flora, and dominant species are R. 
pseudoacacia, Pinus tabulae, Platycladus orientalis, 
Populus davidiana, Caragana microphylla, Lespedeza 
biaolorcz and Zizyphus jujubamill, and so on. The zonal soil 
type is Cinnamon Soil, and Eluvial Soil or Forest Cinnamon 
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Soil can be seen in remote mountains. Soil pH value is in 
range of 5.9-8.5. 

Experimental and measuring method 

Surface runoff fields 

Two fields of surface runoff were separately set up in R. 


pseudoacacia stand and its clearcut area. The natural con¬ 
ditions of the two fields were similar except vegetation 
coverage (Tablel). Intercepting ditch and triangular meas¬ 
uring weir were built at one end of the two fields. The de¬ 
signed discharge of the triangular measuring weir was 
0.212 m 3 -s'\ 


Tablel. Natural conditions of the surface runoff fields 


Land type 

Area 

/hm 2 

Slope 

Direction 

Soil thickness 

/cm 

Matrixes 

Vegetations 

Robinia. 
pseudoaca¬ 
cia stand 

1.00 

25°-30° 

East 

Cinnamon 

40-50 

Calcareousness 

8-10 years old R. pseudoacacia, with coverage of 
80%~90%. Shrub and herbage include Caragana micro- 
phylla and Setaria lutescens, with a coverage of 15% 

Clearcut 

area 

1.04 

20°-28° 

East 

Cinnamon 

30-40 

Calcareousness 

Clearcut was made 2 years ago. Shrub and herbage include 
Setaria lutescen, and Artemisia sp., with a coverage of 30% 


Runoff and sediment 

The runoff volume was measured by sediment pool, 
automatic water-stage recorder, and triangular measuring 
weir. Sediment volume was the sum of dry weights of sus¬ 
pended load and fractional load, while dry weight of sus¬ 
pended load was the weight of sand being filtrated from the 
runoff sample, and that of fractional load was the weight of 
sand depositing in sediment pool (Khai 1997). 

Soil moisture 

The soil moisture contents of 0-10 cm, 10-20 cm and 
20-30 cm layer were separately measured by using the 
oven dry method at the spots on upper, middle and down 
slopes on 1st, 11th and 21st day of each month. At the 
same time, the dynamic laws of soil moisture were ana¬ 
lyzed by the harmonic analysis method. The basic principle 
of harmonic analysis method is to use the known periodic 
function (or time series) and Fourier series method to gain 
the discontinuous spectrum, and then to calculate the 
harmonic factor of key period. According to Fourier theorem, 
under certain conditions, the arbitrary function [f(t)] can be 
expressed by the sum of arbitrary sine (or cosine) trigo¬ 
nometric function. 

/(O = ~ + Z c * ■ sin( k ■ + ^.<>) 

- < -1 

= 'T' + Z ' cos ( k ' •' + 0 t o ) 

- 4-1 


= — + [a* cos( k • CO • /) + b k sin( k ■ co ■ /)] 

- i .1 

where, a k , b k is Fourier coefficient, k is wave numbers of 
harmonic component (k= 0, 1,2...), C k is amplitude of the 

“k” wave; 6 k 0 is first phase of the “k" wave, Ct) =2 n /t is 

basic frequency, t is basic period, t/k is period of the “k” 
wave; f(t) is the arbitrary function; a 0 is Fourier coefficient 
when k is zero. 

Meanwhile, with the random probability decision method 
of A*Schuster, the significant test of the harmonic period 
was carried out. When the periodic ampli¬ 
tude C : - C „ : > 0 - it was significant, and otherwise it was 

not significant (Huang 1983; Guan 1989). 

Results and analysis 

Contrastive analysis of runoff and sediment 

There existed a significant difference in volumes of runoff 
and sediment between the sites of R. pseudoacacia stand 
and its clearcut area. The runoff volume and sediment 
volume in clearcut area were much bigger than those in R. 
pseudoacacia stand (Table 2), with an increase amount of 
40%-177% for runoff and 180%-400% for sediment. This 
indicated that R. pseudoacacia stand could reduce surface 
runoff and oil erosion and increase infiltration of rainfall (Yu 
1991). 


Table2. Runoff volume and sediment volume in R pseudoacacia stand and its clearcut area 


Number 

Amount of 

rainfall 

/mm 

Intensity of 
rainfall intensity 
/mm-min'' 


Runoff volume /mm 


Sediment volume /t ■ 

km' 2 

Wood 

land 

Clearcut area 

Amount of 

increase /% 

Wood 

land 

Clearcut area 

Amount of 

increase /% 

N 0 , 

54.1 

0.55 

7.033 

13.918 

97.9 

30.01 

150.08 

400.1 

No2 

30.5 

0.22 

3.051 

4.491 

47.2 

27.95 

78.34 

180.3 

Noo 

32.4 

1.10 

4.212 

10.972 

160.5 

14.34 

63.58 

343.4 

No4 

17.8 

0.21 

1.780 

3.115 

75.0 

12.90 

38.57 

199.0 

Nos 

42.4 

0.77 

5.310 

12.930 

143.5 

20.00 

100.00 

400.0 

No6 

38.2 

0.84 

4.966 

13.756 

177.0 

21.95 

105.25 

379.5 
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Reason and hydrologic process of runoff occurrence 

Based on the runoff characteristics in hilly semi-arid area, 
the following hydrologic equation was used to express run¬ 
off volume, 

R s (t)= [ idt- [/„ dt- f e-dt -[ s d -dt-\f-dt (1) 

Jo Jo Jo Jo Jo 


arcut area were significantly different (the variance analysis 
results: F=7.823 > F 0 .oi=7.560). The soil moisture in R. 
pseudoacacia stand was higher than that in clearcut area 
with an increase amount of 2.3% (Table 3), especially in the 
rainy season, it was 2.40%-3.66% higher than that in cle¬ 
arcut area. The results demonstrated that R, pseudoacacia 
stand could intercept and reallocate rainfall runoff and in¬ 
crease soil moisture. 


where, R s (t) is amount of runoff on sloping land, i is rainfall 
intensity, i n is rainfall interception rate, e is evaporation 
rate, S<jis filling lowland rate, f is soil infiltration rate. 

During a rainfall, i n , e and S d can be ignored be¬ 
cause they are small and do not cause runoff directly. 
Thus, the amount of runoff depends on the first item (/) and 
last item (/) of equation (1). By calculating equation (1) we 
can obtain Equation (2): 

dR (t) / dt = a = i - f (2) 

From equation (2), we can find the runoff rate of wood¬ 
land mainly depends upon i and f. Under natural rainfall 
condition, runoff occurs when i >t but no runoff occurs 
when i ^ f. 

Therefore, through hydrograph of surface runoff of typical 
rainfall (Figs.1, 2), we could find that the peak value of 
runoff in R. pseudoacacia stand was decreased by 
1.0-2.5x10' 3 m 3 -s'' compared with that in its clearcut area. 
This demonstrated that R. pseudoacacia stand could de¬ 
crease the peak value and volume of runoff. Meanwhile, the 
occurring time of peak value of runoff in R. pseudoacacia 
stand was 10-20 min later than that in its clearcut area, 
which indicated that R. pseudoacacia stand could store 
water and delay the occurring time of flood peak. 

R.pseudoacacia stand 
Clearcut area 


♦ » -V. 

«♦ \ 

* / \ \ 

/ 

J 7 1 ■ 1 2-1^1_i- 

20:03 20:10 20:17 20:24 20:31 20:39 20:46 
Time /hour 

Fig. 1 Hydrograph of surface runoff of the Not typical rainfall 

Harmonic analysis of soil moisture dynamics 

Zero wave value ( 9 0 -the mean soil moisture content) 

The zero wave values of harmonic analysis of soil mois¬ 
ture dynamics (6 0 ) in R. pseudoacacia stand and its cle- 


10 


2 



Time/hour 

Fig. 2 Hydrograph of surface runoff of the N 02 typical rainfall 

Periodic test of harmonic analysis of soil moisture dynamics 

The periodic test S value of harmonic analysis of soil 
moisture dynamics showed that the first waves of soil 
moisture in R. pseudoacacia stand and its clearcut area 
were significant ( <?=0.1), and the second and third waves 
were not significant (Table 4). In harmonic analysis, the first 
wave, the second wave, and third wave are separately 
called the year wave, half a year wave and seasonal wave, 
and they form the key period of time series variation. 
Hence, the results of harmonic analysis also indicated that 
the soil moisture both in R. pseudoacacia stand and its 
clearcut area could be divided into dry season and humid 
season and varied periodically with annual rainfall precipi¬ 
tation, i.e., the soil moisture was lower between the end of 
April to the end of June, began to increase from the early of 
July to mid-September, and decline from the end of Sep¬ 
tember. 

Conclusions and suggestions 

The results of this study showed that pseudoacacia stand. 
The runoff volume and sediment volume in clearcut area 
were much bigger than those in R. pseudoacacia stand 
(Table 2), with an increase amount of 40%-177% for runoff 
and 180%-400% for sediment These indicated that R. 
pseudoacacia stand could decrease surface runoff and soil 
erosion, and increase infiltration of rainfall. 

Through hydrograph of surface runoff of typical rainfall, 
we found that R. pseudoacacia stand could decrease the 
peak value and the runoff volume, store water and delay 
the occurring time of flood peak. 

By the harmonic analysis of soil moisture dynamics, we 
found the soil moisture in R. pseudoacacia stand and its 
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clearcut area were both divided into dry season and humid 
season and varied periodically with annual rainfall precipi¬ 
tation. The soil moisture was higher between the early of 
July to mid-September, and then began to decline from the 
end of September. The soil moisture was lower from the 
end of April to the end of June. The annual mean soil 


moisture in R. pseudoacacia stand was higher than that in 
its clearcut area, in the rainy season, more than 
2.40%-3.66%. So the results indicated that R. pseudoaca¬ 
cia stand could intercept and reallocate rainfall runoff, and 
increase soil moisture. 


Table 3. The zero wave values of harmonic analysis of soil moisture dynamics (0 0 ) 


Land 

Locality 

type 


Soil moisture dynamics ($ ) 


0-10 

Different soil layer /cm 

10-20 

20-30 

Average 

_ . Slope base 

11.6 

10.3 

9.7 

10.5 

R. psudoacacia 





Mid slope 

11.9 

10.6 

10.4 

10.9 

Crest slope 

12.5 

11.8 

11.3 

12.1 

Average 

12.0 

10.9 

10.5 

11.3 

Slope base 

9.2 

8.9 

8.4 

8.8 

Clearcut area Mid slope 

9.2 

8.6 

8.4 

8.8 

Crest slope 

9.8 

9.3 

9.4 

9.5 

Average 

9.4 

8.9 

6.7 

9.0 


Table 4. The periodic test S value of harmonic analysis of soil moisture dynamics 


Land 

Locality 

S value of different locality 


S value of different land type 


type 


First wave 

Second wave 

Third wave 

First wave 

Second wave 

Third wave 

ft. psudoacacia 

Slope base 

1.320 

-0.547 

-1.311 




stand 

Mid slope 

1.610 

-0.544 

-1.932 

1.612 

-0.533 

-1.741 

Crest slope 

1.468 

-0.563 

-1.915 




Clearcut 

Slope base 

0.378 

-0.313 

-0.454 




Mid slope 

0.614 

-0.481 

-0.857 

0.642 

-0.518 

-0.760 

area 

0.467 

-0.504 

-0.758 




Crest slope 

When the periodic amplitude S= q 2 

- C j > 0 ■ i’ is significant («=0.1), and otherwise it is not significant. 




According to the study, the benefits of soil and water 
conservation of artificial protection forest of R. pseudoaca¬ 
cia were obvious. Therefore, we should speed up the con¬ 
struction of artificial protection forest of R. pseudoacacia in 
hilly semi-arid area in China. Moreover, because of R. 
pseudoacacia being also a key tree for fuel forest, we 
should pay more attention to its cutting style, cutting cycle, 
and regeneration after removal of old growth. 
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